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A search is made for massive highly ionising particles with lifetimes in excess of 100 ns, with the ATLAS experiment at the Large 
Hadron Collider, using 3.1 ptr 1 of pp collision data taken at -\fs — 7 TeV. The signature of energy loss in the ATLAS inner detector 
and electromagnetic calorimeter is used. No such particles are found and limits on the production cross section for electric charges 
6e < \q\ < \le and masses 200 GeV < m < 1000 GeV are set in the range 1 - 12 pb for different hypotheses on the production 
mechanism. 
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1. Introduction 



The observation of a massive long-lived highly ionising par- 
^ ticle (HIP) possessing a large electric charge \q\ » e, where e 
i is the elementary charge, would represent striking evidence for 
Q-iphysics beyond the Standard Model. Examples of putative par- 



tides which can give rise to HIP signatures include Q-balls 



i 



i"^ i stable micro black-hole remnants [2], magnetic monopoles [3] 
and dyons [4]. Searches for HIPs are made in cosmic rays [5] 
and at colliders fl; recent collider searches were performed 
at LEP 10-01 and the Tevatron Ir9l4l2ll. Cross sections and 
iy-} event topologies associated with HIP production cannot be re- 
■^j- liably predicted due to the fact that the coupling between a HIP 
and the photon is so strong that perturbative calculations are 
not possible. Therefore, search results at colliders are usually 
quoted as cross section limits in a range of charge and mass for 
7—1 given kinematics Jst] . Also, for the same reason, limits obtained 
I at different collision energies or for different types of collisions 
i> cannot be directly compared; rather, they are complementary. 
HIP searches are part of a program of searches at the CERN 
Large Hadron Collider (LHC) which explore the multi-TeV en- 
d ergy regime. Further motivation is provided by the gauge hi- 
erarchy problem, to which proposed solutions typically postu- 
late the existence of hitherto unobserved particles with TeV- 
scale masses. HIPs at the LHC can be sought at the dedicated 
MoEDAL plastic-track experiment lfl3ll or, as in this work, via 
their active detection at a multipurpose detector. 

Due to their assumed large mass (hundreds of GeV), HIPs 
are characterised by their non-relativistic speed. The expected 
large amounts of energy loss per unit length (dE/dx) through 
ionisation (no bremsstrahlung) are mainly due to the high parti- 
cle charge, but also due to the low speed. The ATLAS detector 
is well suited to detect HIPs. A HIP with sufficient kinetic en- 
ergy would leave a track in the inner detector tracking system of 
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ATLAS and lose its energy on its way to and through the elec- 
tromagnetic calorimeter, giving rise to an electron-like signa- 
ture. The presence of a HIP can be inferred from measurements 
of the proportion of high-ionisation hits in the inner detector. In 
addition, assuming isolation, the lateral extent of the energy de- 
position in the calorimeter is a sensitive discriminant between 
HIPs and Standard Model particles. 

The ranges of HIP charge, mass and lifetime for which unam- 
biguous conclusions can be drawn are determined by the chosen 
trigger and event selections. The choice of an electromagnetic 
trigger limits the phase space to HIPs which stop in the elec- 
tromagnetic calorimeter of ATLAS. The search is optimised 
for data collected at relatively low instantaneous luminosities 
(up to 10 31 crrr 2 s~'), for which a low (10 GeV) trigger trans- 
verse energy threshold is available. In the barrel region of the 
calorimeter, this gives access to energy depositions correspond- 
ing to HIPs with electric charges down to 6e. Standard electron 
reconstruction algorithms are used, which implies that tracks 
which bend like electrically charged particles are sought. Par- 
ticles with magnetic charge, or electric charge above 17e, are 
not addressed here due to the bending along the beam axis in 
the case of a monopole, and due to effects from delta electrons 
and electron recombination in the active detector at the corre- 
sponding values of energy loss (dE/dx > 2 ■ 10 3 MeV/cm). For 
such types of HIPs, more detailed studies are needed to assess 
and minimise the impact of these effects on the selection effi- 
ciency. The 1000 GeV upper bound in mass sensitivity is deter- 
mined by trigger timing constraints, as a significantly heavier 
HIP (with charge 17<? or lower) would be delayed by more than 
12 ns with respect tojS = 1 when it stops in the electromagnetic 
calorimeter (this corresponds to /3 < 0.3), and would thus risk 
being triggered in the next proton bunch crossing. The search 
is sensitive to HIP lifetimes larger than 100 ns since a particle 
which decays much earlier in the calorimeter (even after stop- 
ping) would spoil the signature of a narrow energy deposition. 



Preprint submitted to Physics Letters B 



January 11, 2013 



2. The ATLAS Detector 

The ATLAS detector B14I1 is a multipurpose particle physics 
apparatus with a forward-backward symmetric cylindrical ge- 
ometry and near An coverage in solid angle 1 15]. A thin super- 
conducting solenoid magnet surrounding the inner part of the 
ATLAS detector produces a field of approximately 2 T along 
the beam axis. 

Inner detector (ID) tracking is performed by silicon-based 
detectors and an outer tracker using straw tubes with particle 
identification capabilities based on transition radiation (Tran- 
sition Radiation Tracker, TRT). The TRT is divided into bar- 
rel (covering the pseudorapidity range I77] < 1.0) and endcap 
(0.8 < \rj\ < 2.0) components. A track gives a typical number 
of straw hits of 36. At the front-end electronics of the TRT, dis- 
criminators are used to compare the signal against low and high 
thresholds. While the TRT has two hit threshold levels, there is 
no upper limit to the amount of ionisation in a straw which will 
lead to a signal lfl6ll . guaranteeing that highly ionising parti- 
cles would not escape detection in the TRT. Rather, they would 
produce a large number of high-threshold (HT) hits along their 
trajectories. The amount of ionisation in a straw tube needed for 
a TRT HT hit is roughly equivalent to three times that expected 
from a minimum ionising particle. 

Liquid-argon sampling electromagnetic (EM) calorimeters, 
which comprise accordion-shaped electrodes and lead ab- 
sorbers, surround the ID. The EM calorimeter barrel (I77I < 
1 .475) is used in this search. It is segmented transversely and 
divided in three layers in depth, denoted first, second, and third 
layer, respectively. In front of the accordion calorimeter a thin 
presampler layer is used to correct for fluctuations of energy 
loss. The typical cell granularity (A77 x A(p) of the EM barrel 
is 0.003 x 0.1 in the first layer and 0.025 x 0.025 in the second 
layer. The signal expected for a HIP in the considered charge 
range lies in a region in time and energy where the electronic re- 
sponse in EM calorimeter cells is well understood and does not 
saturate. The robustness of the EM calorimeter energy recon- 
struction has been studied in detail and pulse shape predictions 
are consistent with the measured signals II 1711 . 

The stopping power of a HIP in the ATLAS detector depends 
on its charge, mass and energy, as well as the material budget 
along its path. Details of the latter are given in Ref. Ill 811 in 
terms of number of radiation lengths Xq, as a function of depth 
and pseudorapidity. The integrated radiation length between the 
interaction point and the exit of the TRT is 0.5 Xq at 77 — and 
1.5 Xq at \tj\ = 1.3. The additional amount of material before 
the first layer of the EM calorimeter is 2.0 Xo at 77 = and 3.5 
Xq at \rj\ = 1.3. The thicknesses of the first, second and third 
EM layers are 4.5 X , 16.5 X and 1.5 X at 77 = and 3 X Q , 20 
Xo and 5 Xo at |?7| = 1.3, respectively. 

3. Simulated Event Samples 

Signal events are generated with the Pythia Monte Carlo 
(MC) event generator 1U9I1 according to the fermion pair pro- 
duction process: p + p — > f + f + X. Ref. i20ll is used for 
the parton distributions of the proton. Direct pair production 




ATLAS Simulation 




PYTHIA 

pp -> ffX 
Vi=7 TeV 

m=200 GeV 
m=500 GeV 
m=1000 GeV 



400 500 
E km [GeV] 



Figure 1 : Distributions of pseudorapidity r\ (top) and kinetic energy E#„ (bot- 
tom) at origin for heavy fermions produced with the Drell-Yan process. The 
latter is given with a requirement of \r\\ < 1.35. The distributions for the three 
different masses are normalised to the same number of entries. 



implies that the HIPs are not part of a jet and thus isolated. A 
Drell-Yan-like production mechanism, modified to take into ac- 
count the mass of the HIP 1 21], is used to model the kinematic 
properties of the HIPs. Generated 77 distributions, as well as ki- 
netic energy (Ekin) spectra in the central region (I77I < 1.35), are 
shown in Figure Q] for the three mass points considered in this 
search. 



An ATLAS detector simulation 112211 based on Geant-4 B23H 
is used, where the particle interactions include secondary ion- 
isation by delta electrons in addition to the standard ionisa- 
tion process based on the Bethe-Bloch formula. A correction 
for electron-ion recombination effects in the EM calorimeter 
(Birks' correction) is applied, with typical visible energy frac- 
tions between 0.2 and 0.5 for the signal particles considered. 
Effects of delays are simulated, except for the ability to trigger 
slow-moving particles within the proton bunch crossing time, 
which is considered separately as a systematic uncertainty (see 
Section [5). Samples of approximately 20000 events are pro- 
duced for HIPs with masses of 200, 500 and 1000 GeV. For 
each mass point, HIPs with charges 6e, lOe and 17e are simu- 
lated. 

A data-driven method is used in this work to estimate back- 
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grounds surviving the final selections (see Section I4T2V How- 
ever, in order to demonstrate that the distributions of the rele- 
vant observables are understood, a sample of simulated back- 
ground events is used. The background sample, generated with 
Pythia lfl9ll and labeled "Standard Model", consists mostly of 
QCD events in which the hard subprocess is a strong 2-to-2 
process with a matrix element transverse momentum cut-off of 
15 GeV, but also includes contributions from heavy quark and 
vector boson production. A true transverse energy larger than 
17 GeV in a typical first level trigger tower is also required. 
This sample contains 4 • 10 7 events and corresponds roughly to 
an integrated luminosity of 0.8 pb -1 . 

4. Trigger and Event Selection 

The collected data sample corresponds to an integrated lu- 
minosity of 3.1 + 0.3 pb -1 , using a first level trigger based on 
energy deposits in the calorimeters. At the first level of the trig- 
ger, so-called trigger towers with dimension ArjxA<p = 0.1x0.1 
are defined. In each trigger tower the cells of the electromag- 
netic or hadronic calorimeter are summed. EM clusters with 
fixed size A/7 x Acp = 0.2 x 0.2 are sought and are retained 
if the total transverse energy (£Y) in an adjacent pair of their 
four trigger towers is above 5 GeV. Further electron-like higher 
level trigger requirements are imposed on the candidate, includ- 
ing Ej > 10 GeV, a matching to a track in the ID and a veto 
on hadronic leakage 112411 . The efficiency of this trigger for the 
data under consideration is measured to be (94.0 ± 1.5)% for 
electrons with Ej > 15 GeV and is well described by the sim- 
ulation. The simulation predicts that a highly charged particle 
which stops in the EM barrel would be triggered with a similar 
efficiency or higher. 

Offline electron candidates have cluster sizes of Arj x A<p = 
0.075 x 0.175 in the EM barrel, with a matched track in a win- 
dow of Ar] x A<p = 0.05 X 0.1 amongst reconstructed tracks 
with transverse momentum larger than 0.5 GeV. Identification 
requirements corresponding to "medium" electrons 025]], im- 
plying track and shower shape quality cuts, are applied to the 
candidates. These cuts filter out backgrounds but have a negli- 
gible impact on the signal, for which the cluster width is much 
narrower than for typical electrons. The cluster energy is es- 
timated correcting for the energy deposited outside the active 
calorimeter regions, assuming an EM shower. 

Further offline selections on the cluster transverse energy 
(Ej > 15 GeV) and pseudorapidity (I77I < 1.35) are im- 
posed. The Ej selection guarantees that the trigger efficiency 
is higher than 94% for the objects under study. The restric- 
tion of \tj\ < 1.35 excludes the transition region between the 
EM calorimeter barrel and endcap, reducing the probability for 
backgrounds to fake a narrow energy deposition. 

4.1. Selection Cuts 

A loose selection based on TRT and EM calorimeter infor- 
mation is also imposed on the candidates to ensure that the 
quality of the track and cluster associated to the electron-like 
object is good enough to ensure the robustness of the HIP se- 
lection variables, and to provide a data sample with which to 



estimate the background rate. Only candidates with more than 
10 TRT hits are retained. In addition to the Ej > 15 GeV 
cut for the EM cluster associated with the candidate, a signif- 
icant fraction of the total cluster energy is required to be con- 
tained in six calorimeter cells among the first and second EM 
layers. This is done by requiring the summed energy in the 
three most energetic cells in each of the first and second layers 
to be greater than 2 and 4 GeV, respectively. Following these 
selections, 137503 candidates remain in the data. 

Two sets of observables are used in the final selection. The 
ID-based observable is the fraction, /ht, of TRT hits on the 
track which pass the high threshold. The calorimeter-based dis- 
criminants are the fractions of energies outside of the three most 
energetic cells associated to a selected EM cluster, in the first 
and second EM calorimeter layers: w\ and w%. 

The ff/T distribution for loosely selected candidates is shown 
in Figure [2] The data extend up to far =0.8. The prediction 
of the signal simulation for a HIP of mass 500 GeV and charge 
lOe is also shown. It peaks at fm ~ 1 and has a small tail 
extending into the Standard Model region. 



en 10° 

CD 

1 10 5 k 
° 10 3 k 

<D 
_Q 

E 10 2 
10k 

1 



10 



n I I I r 
ATLAS 

\[s=7 TeV 



• data 3.1 pb" 1 
ES Standard Model MC 

Signal, |q|=10e, 

m=500 GeV 



0.2 



0.4 



0.6 



0.8 



Figure 2: Distribution of the fraction of TRT high-threshold hits for candidates 
satisfying the loose selection. Data (dots) are compared with area-normalised 
signal (\q\ = lOe and m = 500 GeV, dashed line) and Standard Model back- 
ground (shaded area) simulations. The dotted line shows the selection cut value. 

The distributions of w\ and W2 also provide good discrimina- 
tion between signal and background, as shown in Figure [3] For 
a signal, energy is deposited only in the few cells along the par- 
ticle trajectory (as opposed to backgrounds which induce show- 
ers in the EM calorimeter) and the distributions peak around 
zero for both variables. The shapes of the measured distribu- 
tions are well described by the background simulation. A faint 
double-peak structure is visible in data and in background sim- 
ulations for the /ht, W\ and W2 distributions in Figures [2] and 
[3] where the main peak (closest to the signal) corresponds to 
electrons and the secondary peak corresponds to hadrons which 
fake the electron identification signature. 

Finally, the following HIP selection is made: /ht > 0.65, 
w\ < 0.20 and W2 < 0.15. For signal particles, these cuts re- 
ject only candidates in the tails of the distributions, and vary- 
ing them has a minor impact on the efficiency; this feature is 
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probability of a background candidate passing the TRT require- 

and the probability to pass the 



ment is then P 



In 



calorimeter requirements is P w = 



(Ni+N fHT ) 

(ATi+iv,,.) ' leadin 8 t0 an eX " 

pectation of the number of background candidates entering the 
signal region: Nbg = Ni oose P f HT P w - The data sample yields 
N = 0, Ni = 137342, N fliT = 18 and N w = 143, leading 
to P fHT = (1.3 ± 0.3) • 10 4 and P w = (1.0 ± 0.1) • 10~ 3 . The 
expected number of background candidates surviving the selec- 
tion, and thereby the expected number of background events, is 
thus Nb g = 0.019 ± 0.005. The quoted uncertainty is statistical. 
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Figure 3: Distributions of w\ and w% following the loose selection. Data (dots) 
are compared with area-normalised signal (|^| = lOe and m = 500 GeV, dashed 
lines) and Standard Model background (shaded area) simulations. Negative 
values are caused by pedestal fluctuations. Dotted lines show the selection cut 
values. 



Figure 4: Contours of W2 versus fur distributions following loose selection, 
showing the density of entries on a log scale. Data and signal Monte Carlo 
(\q\ = lOe and m = 500 GeV) are shown, and no candidates in the data appear 
near the signal region. The correlation factor between w% and fyr in the data is 
positive (coefficient 0. 15); the same trend is also true for the correlation between 
wi and/z/7- (coefficient 0.18). 



common to all considered charge and mass points. The cut val- 
ues were chosen to yield a very small (<K 1 event) expected 
background (see Section l4~2i > while retaining a high (~ 96%) 
efficiency for the signal. No candidates in data or in simulated 
Standard Model events pass this selection. 

4.2. Data-driven Background Estimation 

A data-driven method is used to quantify the expected back- 
ground yield after the HIP selection. Potential backgrounds 
consist mainly of electrons. For Standard Model candidates, 
the ID and calorimeter observables are correlated in a way that 
further suppresses the backgrounds (see Figure 0. The back- 
ground estimation assumes that /ht is uncorrected with w\ and 
W2 and is thus conservative. 

The yield of particle candidates passing the loose selection 
N/oose = 137503 can be divided into the following: A^o, N\, 
N/ HT , and N w , which represent the number of candidates which 
satisfy both of the selections, neither of the selections, only the 
/ht selection, and only the w\ and W2 selections taken together, 
respectively. Even in the presence of a signal, Ni, N/ HT and 
Af vv . would be dominantly composed of background events. The 



5. Signal Selection Efficiency 

5.7. Efficiencies in Acceptance Kinematic Regions 

The probability to retain a signal event can be factorised in 
two parts: acceptance (probability for a HIP in a region where 
the detector is sensitive) and efficiency (probability for this HIP 
to pass the selection cuts). The acceptance is defined here as 
the probability that at least one signal particle will be in the 
range \rj\ < 1.35 and stop in the second or third layer of the EM 
calorimeter. If this condition is satisfied, the simulation pre- 
dicts a high probability to trigger on, reconstruct and select the 
event. This corresponds to the dark region in Figure |5] which 
shows the predicted selection efficiency mapped as a function 
of the initial HIP pseudorapidity and kinetic energy, in the case 
of \q\ — lOe and m = 500 GeV. Such acceptance kinematic re- 
gions can be parametrised with three values defining three cor- 
ners of a parallelogram. These parameters are summarised in 
Table [TJ For HIPs produced inside such regions, the candidate 
selection efficiency is flat within 10% and takes values between 
0.5 and 0.9 depending on the charge and mass (see Tabled. For 
\q\ — 17e, the main source of inefficiency is the requirement on 
the number of TRT HT hits, which contributes up to 20% signal 
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200 


100 
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500 
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260 
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Table 1: Kinetic energies (in GeV) defining the acceptance kinematic ranges 
for HIPs with the masses and electric charges considered in this search. The 
three columns correspond to the lower left, lower right, and upper left corners 
of parallelograms in the (\r]\, E^in) plane. 



m [GeV] 


\q\ = 6e 


\q\ = lOe 


\q\ = lie 


200 


0.822 ± 0.026 


0.820 ±0.015 


0.484 ±0.012 


500 


0.868 ± 0.021 


0.856 ±0.014 


0.617 ±0.011 


1000 


0.558 ±0.019 


0.858 ±0.012 


0.700 ±0.012 



Table 2: Expected fractions of HIP candidates passing the final selection, as- 
suming they are isolated and produced inside the acceptance regions defined 
by the values in Table [T] Uncertainties due to MC statistics are quoted; other 
systematic uncertainties are discussed in Section|6] 



loss. This is largely due to the presence of track segments from 
delta electrons, which have a non-negligible probability to be 
chosen by the standard electron track matching algorithm. For 
low charges, inefficiencies are dominated by the cluster Ej cut, 
typically accounting for ~ 6% loss. Other contributions, like 
trigger, electron reconstruction, and electron identification, can 
each cause 1 - 6% additional inefficiency. 



5.2. Efficiencies for Drell-Yan Kinematics 

The estimated fractions of signal events where at least one 
candidate passes the final selection, assuming they are produced 
with Drell-Yan kinematics, are shown in Table |3]for the values 
of charge and mass considered in this search. The dominant 
source of loss (70 - 85% loss) is from the kinematic accep- 
tance, i.e., the production of HIPs with I77I > 1.35, as well as 
their stopping before they reach the second layer of the EM 
calorimeter, or after they reach the first layer of the hadronic 
calorimeter. The relative contributions from these various types 
of acceptance loss depend on mass and charge, as well as the 
kinematics of the assumed production model. The Drell-Yan 
production model implies that the fraction of HIPs produced in 
the acceptance region of pseudorapidity \rj\ < 1 .35 is larger with 
increasing mass (see Figure [1}. Also, with the assumed energy 
spectra (bottom plot in Figure Q3, the acceptance is highest for 
intermediate charges (\q\ — lOe), since HIPs with low charges 
tend to punch through the EM calorimeter and HIPs with high 
charges tend to stop before reaching it. 




Figure 5: Probability to pass all selec tion criteria as a function of pseudorapidity 
and kinetic energy at origin, for a HIP with charge lOe and mass 500 GeV. The 
dark region corresponds to the kinetic range where the particle stops in or near 
the second layer of the EM calorimeter barrel and is parametrised with three 
energy values (dashed parallelogram, see Table^. 



m [GeV] 


\q\ = 6e 


\q\ = lOe 


\q\ = Vie 


200 


0.102 ±0.002 


0.175 ±0.003 


0.1 12 ±0.002 


500 


0.150 ±0.003 


0.236 ± 0.003 


0.193 ±0.003 


1000 


0.133 ±0.002 


0.299 ± 0.004 


0.237 ± 0.004 



Table 3: Expected fractions of signal events passing the final selection, assum- 
ing Drell-Yan kinematics. Uncertainties due to MC statistics are quoted; other 
systematic uncertainties are discussed in Section [6] 



6. Systematic Uncertainties 

The major sources of systematic uncertainties affecting the 
efficiency estimation are summarised below. These mainly con- 
cern possible imperfections in the description of HIPs in the 
detector by the simulation. 

• The recombination of electrons and ions in the sampling 
region of the EM calorimeter affects the measured current 
and thus the total visible energy. Recombination effects 
become larger with increasing d£7dx. In the ATLAS simu- 
lation, this is parametrised by Birks' law ||26j. To estimate 
the uncertainty associated with the approximate modeling 
of recombination effects, predictions from the ATLAS im- 
plementation of Birks' correction Il27ll are compared to 
existing data of heavy ions punching through a layer of 
liquid argon 1I28U30I1. In the range 2 • 10 2 MeV/cm < 
d£7dx < 2- 10 3 MeV/cm, which corresponds to typical HIP 
energy losses in the EM calorimeter for the charges and 
masses under consideration, the uncertainty in the sim- 
ulated visible energy fraction is ±15%. This introduces 
between 4% and 23% uncertainty in the signal selection 
efficiency. The impact is largest for charge 6e, for which 
a lower visible energy would be more likely to push the 
candidate below the 15 GeV cluster Ej threshold. 



• The fraction of HIPs which stop in the detector prior to 
reaching the EM calorimeter is affected by the assumed 
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amount of material in the Geant-4 simulation. Vary- 
ing the material density within the assumed uncertainty 
range (± ~ 10% MM), independently in the ID and EM 
calorimeter volumes, leads to a 6% uncertainty in signal 
acceptance. 

• The modeling of inactive or inefficient EM calorimeter re- 
gions in the simulation results in a 2% uncertainty in the 
signal efficiency. 

• Cross-talk effects between EM calorimeter cells affect the 
wi and W2 variables and this may not be accurately de- 
scribed by the simulation for large energy depositions per 
cell. The resulting uncertainty in signal efficiency is 2%. 

• Secondary ionisation by delta electrons affects the track 
reconstruction and the calorimeter energy output. The 
amount of delta electrons in ATLAS detectors as described 
in Geant-4 depends on the cutoff parameter (the radius be- 
yond which delta electrons are considered separate from 
the mother particle). Varying this parameter results in a 
3% uncertainty in the signal efficiency. 

• For clusters delayed by more than 10 ns with respect to the 
expected arrival time of a highly relativistic particle, which 
corresponds to /3 < 0.37, there is a significant chance that 
the event is triggered in the next bunch crossing by the first 
level EM trigger. In most of the mass and charge range 
considered in this search, more than 99% of the particles 
which are energetic enough to reach the EM calorimeter 
and pass the event selection are in the high-efficiency range 
P > 0.4. The only exception is \q\ — be and m - 1000 GeV, 
for which the f3 distribution after selection peaks between 
0.32 and 0.47. The trigger efficiency loss is corrected for, 
resulting in an additional 25% uncertainty for this particu- 
lar case. 

• Uncertainties in the choice of parametrisation for the par- 
ton density functions (pdfs) of the proton have an impact 
on the event kinematics. To test this effect, events were 
generated (see Secti on |3] l with 7 different pdfs from vari- 
ous sources |2(1 32-3511 Assuming that acceptance vari- 
ations due to the choice of pdf are Gaussian, the resulting 
relative uncertainty in the acceptance is 3%. 

• The relative uncertainty in efficiency due to MC statistics 
is of the order of 2%. 

Other effects, like event pile-up and electron pick-up by pos- 
itively charged particles, have been investigated and found to 
be negligible. Efficiency systematics are dominated by Birks' 
correction. The relative uncertainties in the signal selection ef- 
ficiencies (Tables [2] and [3]), obtained by adding all effects in 
quadrature, are shown in Table |4] 

The systematic uncertainty in the absolute integrated lumi- 
nosity is 11% Q. 



m [GeV] 


\q\ = 6e 


\q\ = lOe 


\q\ = lie 


200 


25% 


11% 


9% 


500 


17% 


10% 


9% 


1000 


28% 


10% 


9% 



Table 4: Relative systematic uncertainties in efficiency, combining in quadra- 
ture all the effects described in the text. 



to [GeV] 


\q\ = 6e 


\q\ = lOe 


\q\ = lie 


200 


1.4 


1.2 


2.1 


500 


1.2 


1.2 


1.6 


1000 


2.2 


1.2 


1.5 



Table 5: Inclusive HfP cross section upper limits (in pb) at 95% confidence level 
for isolated long-lived massive particles with high electric charges produced in 
regions of pseudorapidity and kinetic energy as defined in Table[T] Efficiencies 
in Table ff] and uncertainties in Tabled were used in the cross section limit 
calculation. 



7. Upper Limit on the Cross Section 

A very low (<K 1 event) background yield is expected and 
no events are observed to pass the selection. Knowing the in- 
tegrated luminosity (3.1 pb~') and the selection efficiency for 
various model assumptions (Tables|2]and|3]l, cross section limits 
are obtained. This is done using a Bayesian statistical approach 
with a uniform prior for the signal and the standard assumption 
that the uncertainties in integrated luminosity (11%) and effi- 
ciency (Table HI are Gaussian and independent. The limits are 
presented in Table [5] (for a particle produced in the acceptance 
kinematic region defined by Table [TJ and in Table [6] (assuming 
Drell-Yan kinematics). 

These limits can be approximately interpolated to intermedi- 
ate values of mass and charge. Also, the limits quoted in Ta- 
ble [5] can be used to extract cross section limits for any given 
model of kinematics by correcting for the acceptance (fraction 
of events with at least one generated HIP in the ranges defined 
by TableQ]): such a procedure yields conservative limits thanks 
to the fact that candidates beyond the sharp edges of the accep- 
tance regions defined in TableQ]can also be accepted. 



8. Summary 

A search has been made for HIPs with lifetimes in excess 
of 100 ns produced in the ATLAS detector at the LHC using 
3.1 pb~' of pp collisions at -\fs = 7 TeV The signature of 



to [GeV] 


\q\ = 6e 


\q\ = lOe 


\q\ = 17e 


200 


11.5 


5.9 


9.1 


500 


7.2 


4.3 


5.3 


1000 


9.3 


3.4 


4.3 



Table 6: Pair production cross section upper limits (in pb) at 95% confidence 
level for long-lived massive particles with high electric charges, assuming a 
Drell-Yan mechanism. Efficiencies in Tableland uncertainties in Table^were 
used in the cross section limit calculation. 
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high ionisation in an inner detector track matched to a narrow 
calorimeter cluster has been used. Upper cross section lim- 
its between 1.2 pb and 11.5 pb have been extracted for HIPs 
with electric charges between 6e and 17e and masses between 
200 GeV and 1000 GeV, under two kinematics assumptions: a 
generic isolated HIP in a fiducial range of pseudorapidity and 
kinetic energy, or a Drell-Yan fermion pair production mech- 
anism. HIP mass ranges above 800 GeV lllE are probed for 
the first time at a particle collider. These limits are the first 
constraints obtained on long-lived highly charged particle pro- 
duction at LHC collision energies. 
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